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ABSTRACT 

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within 
the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint 
Herschel sources near the confusion noise. The method is based on the use of the positions of 5p;Yzer/MIPS 24 /jm sources as priors. Our catalogs 
are able to reliably (5cr) recover galaxies with fluxes above 6 and 10 mjy in the PACS 100 and 160 yum channels, respectively, and 12 to 18 mjy 
in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the 
Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or 
dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric 
redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments. 
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1. Introduction 

Based on the IRAS, ISO, and Spitzer missions, we know 
that luminous infrared galaxies (LIRGs) experienced signif- 
icant evolution in the past 8 Gyr, roughly decreasing their 
typical infrared luminosit y rL(IR)1 by an o r der of magni- 
tude fr om z~ l to z=0 (Floresetal.' '1999'; 'Chary & Elbaz? 
20011: iPerez- Gonzalez et al. 2005; Le Floc'h et al. 2005; 
Rodighiero et all l2010t iMaenelli et all l2009h . By z~l, they 
dominated the star formation activity of the Universe, being 
responsible for more than 50% of the cosmic SFR density up 
to at least z~3, and playing an important role in the formation 
of massive g alaxies at z>1.5, which exhib it high star formation 
efficiencies ( iPerez-Gonzalez et all 120081: iReddv et al.1 l2008t 
[Damen et al.' l2009| ), strong dust ob scured nuclear activity 
(Daddietal. 2007 riDonlev et al.ll2008h . or both. Beyond z~l, 
sub-mm surveys have also detected a population of LIRGs 
and (mostly) ULIRGs, identified as very active stages in the 



formation of massive galaxies (jS mail et al. 1997 ;_Hughes et al] 
ll998HChaDman et al.ll2003[ 120051: , Pone et al.ll2006i) . 

Although very useful for measuring the obscuration of op- 
tical light by interstellar dust (very abundant in the most ex- 
treme starbursts), the analysis of IR data to understand the for- 
mation and evolution of galaxies in the Universe has been clas- 



sically hampered mainly by two issues. First, the sensitivity 
of IR telescopes has increased significantly from one mission 
to the next (and also their angular resolution, although at a 
lower rate), but there are still problems correctly identifying all 
sources detected in IR surveys, most noticeably at longer wave- 
lengths where confusion remains severe (e.g., at 70/160 //m in 
the Spitzer/MlPS surveys and in the SCUBA or LABOCA sub- 
mm data). Moreover, finding counterparts and getting reliable 
redshifts from spectroscopy for high redshift (U)LIRGs is chal- 
lenging since a good fra ction of them are extremely faint at 
optical/NIR wavel engths dlvison e t al. 2002; Wilson et al. 200|; 
Fra ver et al.ll2004l) . A second problem for the study of 0<z<4 
(U)LIRGs is the availability of only one photometric data point 
(MIPS 24 jum) for most of them. The fits to dust emis sion mod- 
els to calculate L(IR), and from this derive a SFR (Kennicutj 
1998), need a considerable extrapolation, and t hus the SFR is 
subje ct to significant systematic un certainties dPapovich et al.l 
.20061: iDaddiet aI.ll2007HBarro et alJUoiQ) . In addition, the lim- 
ited FIR data already available for galaxies at z>2 indicate 
that they have different dust properties from those observed for 
nearby (U)LIRGs (Chapman et al. 2005; Papovich et al. 200l 
iRigbv et alJl2008l:lRieke et alJl2009h~ 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



The ESA Herschel Space Observatory dPilbratt et al.ll2010h 
alleviates these two issues. Thanks to its sensitivity and angu- 
lar resolution at 100</l<500yum, it enables us to robustly char- 
acterize high-z galaxies. In addition to the populations of IR 
galaxies newly discovered by Herschel, it will supply up to 
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six photometric data points in the IR range for 0<z<4 galax- 
ies detected by MIPS at 24 jim. These data can be used to 
constrain the fits to the dust emission templates. Even more 
novel and relevant is the combination of Herschel observations 
with the MIPS data (and also IRAC) to diminish the effects 
of source confusion in IR surveys. The new Herschel data in- 
deed allow an easier and more reliable identification of counter- 
parts from one band to the adjacent one as we move to shorter 
wavelengths. In addition, the Herschel bands effectively fill 
the gap between MIPS and (sub-)mm surveys (with SCUBA2, 
LABOCA, or AzTEC). By combining all these datasets, we can 
now use the several IR photometric data points to estimate a 
photometric redshift based on dust emission alone, and com- 
pare it with the estimations based on UV-to-NIR data. This can 
be regarded as an extension of the radio-to-IR photo-z method 
(Carilli & Yun 1999; Dunne et al. 2000; Renaaraian & Takeuchi 
r2001UBlain et al.ll2003l;IAretxaga et al . 2003;'Daddi et al. 2009), 
benefiting from a finer sampling of the emission from dust at 
different temperatures. The procedure can be used to improve 
the cross-correlation with optical/NIR galaxy samples, or even 
to obtain the only redshift measurement possible for very faint 
(high-z) IR-bright optically-faint sources. These sources may be 
numerous in the Herschel surveys, because of its impressive 
sensitivity at wavelengths probing the peak of the dust emis- 
sion, which can be brighter than the optical emission by orders 
of magnitude. 

In this Letter, we perform a preliminary analysis of the 
Herschel Lensing Survey (HLS, PI Egami) data taken for the 
Bullet cluster to show the utility of position priors and photo- 
metric redshifts obtained from UV-to-NIR and IR-to-mm data to 
improve the identification of Herschel and sub-mm sources with 
those detected at shorter wavelengths. A correct identification of 
IR sources will be extremely relevant to all Herschel cosmolog- 
ical surveys, but is of critical importance for clusters, because of 
their crowded nature and the possible presence of high-z lensed 
and distorted sources with very faint optical counterparts that 
may even only be detectable in the FIR/mm range. 

2. Data and catalogs 

The reduction of the Herschel SPIRE (Griffin et al.] 1201 Oh and 
PACS (^oglitsch et al. 2010) data of the B i illet cl uster (a~06:58 
5 — 55:57) is described in lEgami et alj (1201 Ol) and iRex et all 
( I2010.) . jointly with optical-to-mm ancillary data used to esti- 
mate photometric redshifts, stellar masses, and SFRs. Here we 
describe our cataloging procedure for the five Herschel bands: 
PACS 100 and 160 jum, and SPIRE 250, 350, and 500 yum. 

Given the remarkable depth of the Spitzer MIPS 24 /vm im- 
ages, PACS catalogs have been compiled with a position prior 
technique, in volving a list of MIPS sou rces and a PSF fitting 
analysis (see iPerez-Gonzalez et al. 2008). To improve the reli- 
ability of the method, we first aligned all Herschel images to 
the W CS of the M IPS 24 /^m map using the wcstools soft- 
ware dMinkl Il999l) . obtaining matching uncertainties smaller 
than one pixel for each Herschel band. The MIPS 24 jum cat- 
alog in the Bullet cluster field has a 5cr detection level of 
55o-[24]=85+30 jjiya and a density of 5.3 sources/arcmin^ 
above that threshold (7.2 sources/arcmin^ with S [24]>50 /^Jy). 

The similar PSF sizes of the MIPS 24 //m and the PACS 
100 yum bands (~6" and ~8", respectively) make the identifica- 
tion and extraction of sources straightforward for this band. In 



Table 1. Properties of the HLS Bullet cluster catalogs. 





PACS 100 


PACS 160 


SPIRE250 


SPIRE350 


SPIRE500 


A" 


28.5 


29.7 


434 


438 


416 


S(^)'' 


5.5 (3.0) 


10(6) 


12(7) 


17(10) 


18(11) 


D"^ 


2.0 (5.0) 


2.0 (4.2) 


1.1(1.8) 


0.5(1.2) 


0.2 (0.6) 


M'' 


12% (24) 


19% (24) 


48% (24) 


42% (250) 


48% (350) 



' Calculated from sever al artificial apertures formed with random sky 
pixels; see Appendix A in lPerez-Gonzalez et all l l2008h . 



'' Surveyed area in arcmin^ . 

*" 5cr (3(r) flux density detection levels' in mjy. 

'^ Source densities for each significance level in sources/arcmin^ . 

'' Fraction of merged sources and wavelength (in /jm) used to obtain 

position priors, 
addition, we also considered that some sources may be too faint 
for MIPS yet still be detectable by PACS. To account for this 
population, we detected sources directly in the PACS data using 
Sextractor (Bertin & Arnouts 1996), obtaining a sample of 43 
sources above the 5cr level (55o-[100]=5.5±0.7 mJy). The cross- 
correlation of this sample with the MIPS catalog (using a search 
radius of 3") revealed that all but one of the robust PACS detec- 
tions are detected at 24 jiva, and 5 of them had 2 MIPS counter- 
parts within the search radius. The only non-detection at 24 //m 
is, in fact, a group of three very faint 100 //m sources within 
10", which are detected as one single source by Sextractor in the 
PACS image, but are individual sources at 24 //m. The 5 multiple 
detections were sources for which the slightly higher resolution 
and greater depth of MIPS allowed the deblending of close IR 
emitters. This cross-correlation also found that the WCS align- 
ment between both images has an rms smaller than 1". 

For the 160 jum image, it is expected that more 24 jum 
sources are merged to form single sources. For this reason, be- 
fore measuring photometry, we removed sources from the main 
catalog that could not be separated at the resolution of the PACS 
160 yum band, i.e., sources closer than 6", which is half of the 
PSF FWHM, keeping only the brightest source at 24 yum for 
each group of merged sources. The master catalog combined the 
MIPS source list with the direct detections in the 160 jum image 
above the 5cr level (5 50- [160] =10+1 mJy), all the latter being 
actually detected in the MIPS data. Merged sources account for 
almost 20% of the master catalog, with 85% (15%) of this sub- 
sample having 2 (3) objects within 6". Just considering sources 
directly detected at 160;um, the fraction of them presenting mul- 
tiple identifications in the MIPS catalog is higher than 40%. The 
WCS accuracy between the MIPS and PACS 160 images is infe- 
rior to 1.5". 

The final purged list of sources was used to fit PSFs at the 
given positions using the daophot package in IRAF, allowing 
for one pixel centering offsets. The results of this PSF fitting 
algorithm are flux densities in a given aperture (we found the 
optimum values to be 5" and 7.5" for the PACS green and red 
channels, respectively), to which an aperture correction must be 
applied. The values we applied are 1.84+0.08 and 1.93+0.08 for 
100 and 160jum, respectively. We note that these aperture correc- 
tions are larger than those published by the PACS Team, but we 
found that their PSF matched ours (constructed with the same 
Bullet cluster data) for large radii, but was brighter for similar 
and smaller radii than that selected for the PSF fitting method. 

The SPIRE catalogs were also compiled with a method based 
on priors and PSF fitting. The photometry was carried out in cir- 
cular apertures (of radii 12", 18", and 22" for 250, 350, and 
500 //m, respectively), applying a calibration based on the beam 
sizes (assuming they are Gaussians). For the 250 jum channel, 
we combined the list of sources detected by MIPS at 24 /im and 
the direct detection in the SPIRE data using Sextractor and a 5cr 
threshold (55^[250]=12±2 mJy). Other objects within 9" (haff 
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Fig. 1. Postage stamps of two high-z sources in the Bullet cluster field detected by Herschel: the southern tip of the arc 
MIPS06:58:3 1.1 -55:56:39.2 {left), and the LABOCA source MIPS06:58:45. 3-55:58:46.5 (right). The size of the images is 2'x2', 
except for the HST images, whose size is 40"x40"(shown in the other stamps with a square). Circles (3" radius) pinpoint sources 
detected at 24 /im. Cyan circles show the actual >5cr detections in e ach Herschel band. Gr een circles in the IRAC RGB images 
mark z>2 galaxies in each field, including the three lensed sources in iGonzalez et al.l (l2009l to the SE or the arc) and a LABOCA 
source from Johansson et aLj (.201(1 to the W of the arc). 



of the PSF FWHM) of a given candidate source were purged, 
only the position of the brightest galaxy at 24 yum being kept. 
For the region covered simultaneously by MIPS and the SPIRE 
250 fiva channel, only 2 SPIRE sources out of 261 were not de- 
tected by MIPS, but both of them turned to be F(24)~60 yuJy de- 
tections at the 3cr level. Around half of the SPIRE sources have 
2 or more MIPS counterparts within 9". The WCS alignment 
between the MIPS and SPIRE250 images has an rms of ~4". 

The catalog for the SPIRE 350 ;um channel was constructed 
from the hst of detections at 250 yum. Sources within 14" were 
merged into the same object (about 40% of the total catalog, 
virtually all of them being a doublet), only the position of the 
brightest one at 250 /im being kept. A direct detection in the 
350 yum image (5'5o-[350] = 17+3 mJy) found that 5% of sources 
are not directly linked to a 250 //m emitter, but virtually all of 
these are a combination of several 250 /vm sources that causes 
a displaced 350 /im centroid. Only one isolated high-confidence 
350 /vm source was not detected in the 250 /im catalog, but a 
reliable flux could be recovered in the blue channel once the 
350 //m detection revealed its position. The WCS alignment be- 
tween the MIPS and SPIRE350 images has an rms of ~6". 

For the 500 //m channel (55(r[500] = 18±4 mJy), we started 
from the 350 /^m catalog, and we did not find any object without 
an identification in the bluer bands using a 18" search radius 
(approximately half of the PSF FWHM). Around 50% of the 
SPIRE 500 /im sources were doublets in the 350 //m band. The 
WCS alignment between the MIPS and SPIRE500 images has 
an rms of ~9". 

The prior-based photometric technique was compared with 
the direct blind detection in the Herschel images. In PACS, a 
direct detection down to 5cr included more than 40% of sources 
with no counterpart in the MIPS map, implying that these are 
probably spurious detections. The fluxes of sources extracted 



with both the blind and prior-based detections are virtually iden- 
tical. In SPIRE, the prior-based procedure is able to recover 
fluxes for 30% more sources than the direct detection. For the 
sources in common, the average difference in the photometry is 
2% (the blind method provides brighter fluxes), with a 6% scat- 
ter. 

Table [1] indicates the surveyed areas, detection thresholds, 
the corresponding source densities, and fraction of merged 
sources for our catalogs of the PACS and SPIRE HLS data. An 
important (and potentially dangerous) step in our photometric 
procedure based on position priors is the blending of nearby 
non-resolved sources, which provides only the position of the 
brightest galaxy in the parent catalog to the PSF fitting algo- 
rithm. We note that the MIPS 24 //m flux is dominated by the 
emission of warm dust and/or PAHs while the SPIRE bands are 
dominated by the emission of cold dust. Thus, it is probable that 
bright MIR emitters are not the dominant sources at FIR wave- 
lengths (see Figure [TJ. In the following section, we illustrate a 
method to help in the robust identification and/or deblending of 
the Herschel sources. 

3. Identification of h\Q\\-z Herschel sources 

Figure [1] shows thumbnails in several bands of two interest- 
ing sources detected by Herschel in the Bullet cluste r field: 
the so uthern tip of the arc at z=3.24 described in Mehle rtet al.l 
(l200lh. and a sub-mm ga laxy detected by LABOCA (source #10 



Johansson et al■ll2blO^ . With these two examples, we demon- 



strate that the use of IRAC and MIPS data helps to correctly 
identify the Herschel sources with optical/NIR counterparts. In 
the case of the arc, the photometric method based on the use of 
the MIPS 24 fiva catalog as a prior helped in the detection and 
measurement processes for the PACS 160 //m and SPIRE bands. 
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where the source is very faint and almost indistinguishable from 
the background. In the PACS red band, our method allows a de- 
blending of the source in the two knots detected by MIPS and 
linked to the arc. In the case of the LABOCA source, although 
it is barely detected in the optical, the position of the source can 
be followed as we move to redder wavelengths and larger PSFs, 
allowing a more reliable identification in spite of the presence 
of a close neighbor to the west. The brightest source at 24 yum is 
indeed already brighter than its companion at IRAC wavelengths 
(most noticeably, at 8.0 jum), and this nearby source disappears 
in PACS. This illustrates the benefits of the use of Spitzer data 
as priors to obtain the most probable identifications of Herschel 
and (sub-)mm emitters at shorter wavelengths. 

Figure|2]shows the optical-to-mm SEDs of the two examples 
of high-z galaxies. Photometric redshifts for all 24 fiva sources 
in the Bullet cluster field were estimated from the UV-to-MIR 
data (including IRAC fluxes) usi ng the technique described in 
iPerez-Gonzalez et al.1 (l2005l I2008D . Briefly, the SEDs compiled 
from aperture-matched photometry are fitted with ste llar popu- 
lation models and AGN templates from'P olletta et alj ( |2007|) . A 
minimization algorithm is used to search for the optimal red- 
shifted template fitting the data, and the most probable photo- 
metric redshift is calc ulated by integratin g the redshift probabil- 
ity function (see also iBarro et alJl2010l) . Photometric redshifts 
were also estimated from the IR data alone by comparing the 
Spitzer, Her schel and (sub-)mm fluxes with dust emission tem- 
plates from IChary & Elba3 (l200lh . iDale & HeloiH (l2002h . and 
iRieke et aP (l2009l) . Using these two fits, we estimated relevant 
parameters such as the stellar mass and the SFR. The Rie ke et al] 
(1200 9) template most closely fitting the data for the arc corre- 
sponds to a local L(IR)=10"^ L© galaxy. This value is consistent 
with the average for the arc based on the three template libraries, 
L(IR) =10'^"^'''^"""^ Lq, on ce a magnification factor of ~21 is ap- 
plied dBradac et al.l 120061) . The LABOCA source qualifies as a 
hyper-LIRG based on the best fits to all template libraries (if it 
is not lensed), and the best-fit template from Rieke et al. is that 
of a local galaxy with lO'^^^ L©, although our source seems to 
have more prominent PAH emission. 

Figure|2]illustrates the power of using independent estimates 
of the photometric redshift based on U V-to-NIR data and IR data 
alone to validate the individual values (see the photo-z proba- 
bility distributions). The procedure is also useful for achieving 
more reliable identifications of the IR emitters. For example, for 
the LABOCA source at z~2.5, we calculated photometric red- 
shifts based on UV-to-NIR data for the three close neighbors de- 
tected by MIPS within ~10" of the central source. The galaxy 
to the north lies at z~0.3 (a cluster member), that to the south 
at z~1.9, and the closest companion has a photo-z compatible 
with the central source (z~2.3+0.3), thus implying that the two 
sources are a possible interacting pair. 

The general validity of our method for the whole sample of 
Herschel sources will be tested thoroughly in forthcoming pa- 
pers by evaluating the quality of our photo-z's. From a prelim- 
inary comparison of the optical and FIR photo-z's for sources 
detected by MIPS, PACS, SPIRE, and a (sub)-mm instrument, it 
is encouraging that we measure an average A(z)/(1 + z)-r]-0.09 
with 10% of catastrophic outliers (7;> 0.2), comparable to th e typ- 
ical goodness of IR-based photo-z's dAretxaga et alj|2003l) . 
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Fig. 2. SEDs of the z=3.24 LIRG MIPS06:58:3 1.1 -55:56:39.2 
(magnification factor ~21) and the z~2.5 HyLIRG 
MIPS06:58:45.3-55:58:46.5, detected in the MIR/FIR by 
Spitzer and Herschel, and the (sub-)mm by LABOCA and 
AzTEC. The blue line shows the best-fit stellar population 
model for the photometry up to 8 fj.m, which provides an 
estimate of the photometric redshift and the stellar mass. The 
cyan, green, and red lines show the best-fit dust emission model 
for the MIR- t o-mm phot ometry (from Ch arv & Elbaz 200ll 
Dal e & Heloul I2002L and lRiekeetal.1 l2009l) . which gives an 
independent photometric redshift estimation and a SFR value 
(given in the plots with no correction for magnification factor). 
The inset plots depict the photometric redshift normalized 
probability distributions for the optical/NIR and the MIR/mm 
estimations. 
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